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Forward jet and multijet production has been measured at low Bjorken x at HERA. The 
measured cross sections and correlations were compared to predictions from DGLAP-based 
fixed-order calculations. Further comparisons were made to DGLAP-based and CCFM-based 
leading-order Monte Carlo predictions, as well as to Colour-Dipole model predictions. For the 
majority of the phase space covered in the HERA kinematic region, fixed-order calculations 
describe the data well, while the leading-order models provide an inconsistent description of 
the data. 



1 Introduction 

Jet production in DIS is an ideal environment for investigating different approaches to parton 
dynamics at low Bjorken-x, xbj- An understanding of this regime is of particular relevance 
in view of the startup of the LHC, where many of the Standard Model processes such as the 
production of electroweak gauge bosons or the Higgs particle involve the collision of partons 
with a low fraction of the proton momentum. 

In the usual collinear QCD factorisation approach, the cross sections are obtained as the con- 
volution of perturbative matrix elements and parton densities evolved according to the DGLAP 
evolution equations. In these equations, all orders proportional to as In and terms with 
double logarithms InQ^ • Inl/x, where x is the fraction of the proton momentum carried by 
a parton, which is equal to xbj in the quark-parton model, are resummed. In the DGLAP 
approach, the parton participating in the hard scattering is the result of a partonic cascade 
ordered in transverse momentum, pT- The partonic cascade starts from a low-pT and high-x 
parton from the incoming proton and ends up, after consecutive branching, in the high-pj- and 
low-x parton entering in the hard scattering. At low xbj, where the phase space for parton 
emissions increases, terms proportional to a^lnl/x may become large and spoil the accuracy 
of the DGLAP approach. In this phase-space region, a better description may come from the 
BFKL approach, which resums terms proportional to Inl/x, and the CCFM approach, which 
uses unintegrated gluon densities in an all-loop non-Sudakov resummation. 

Parton evolution schemes at low xbj were studied at HERA by measuring forward jet pro- 
duction and correlations in jet angles and transverse momentum. An excess of forward jets 
compared to DGLAP-based predictions and jets produced in the hard scatter that are not 
strongly correlated in transverse momentum may indicate the breakdown of DGLAP dynamics. 
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2 Dijet Azimuthal Correlations 



Dijet azimuthal correlations were investigated by 
the HI CollaboratioiP by measuring the cross- 
sections d^cj/fixBjdAc/)*, where A(^* is the azimuthal 
separation in the hadronic centre-of-mass (HCM) 
frame between the two selected jets closest to the 
scattered electron in pseudorapidity, rj. The mea- 
surements of A(/)* are reasonably well-described by 
NLOjetEI calculations at 0{a^), albeit within large 
theoretical uncertainties. To reduce the theoretical 
uncertainties, the measurements were normalised to 
the visible cross section for A0* < 170°. With a 
reduced theoretical uncertainty, the calculations are 
shown to predict a narrower A(p* spectrum than is 
measured, especially at very low shown in 

Fig. [TJ The measurements were also compared to 
predictions from two RAPGAplSI (DGLAP) samples, 
with one sample using only direct photons, the other 
using both direct and resolved photons; Lepto^ 
(CDM); and CascadeEI (CCFM). Ah models fail 
to describe A(j)* over the entire range in xbj covered. 
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Figure 1: Double-differential normalised (see 
text) cross sections as a function of Acp* as mea- 
sured by HI compared to NLOjet calculations 
for 0{a^,) and 0{al). 



3 Multijet Correlations 



The sensitiviy of parton evolution to the topology 
of the jet system was studied by the ZEUS collabor- 
tion^. Multi-differential cross sections as functions 
of the jet correlations in transverse momenta, az- 
imuthal angles, and pseudorapidity have been mea- 
sured for dijet and trijet production in the HCM 
frame. DGLAP-based calculations from NLOjet 
at ©(a^) and 0(af ) were compared to the measure- 
ments. The NLOjet calculations at 0(0^) do not 
describe the correlations in transverse momenta and 
azimuthal angle for dijet events; however with in- 
clusion of higher-order terms, the NLOjet calcula- 
tions at 0{a^) describe the dijet data over the entire 
range in xbj covered. The importance of higher- 
order terms at low xbj is seen especially when mea- 
the double-differential cross sections in 
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Figure 2: ZEUS dijet cross sections in xb; and 



as a function of XBj with A(p^^^'^,j < 120' 



^^HCM azimuthal separation of the two jets 

with the highest transverse energy. At low xbj, the 
NLOjet calculations at 0(a|) are up to about one 
order of magnitude larger than the 0{a'g) calcula- 
tions and are consistent with the data, as seen presented in Fig. [2j The NLOjet calculations 
at 0{a^g) also provide a reasonable description of the trijet measurements, with the description 
improving somewhat at higher xbj- 



4 Forward Jet Production 



To examine the sensitivity of parton evolution to forward jet production, the ZEUS cohaboration 
has studied jet production in an extended pseudorapidity range of ?7^^^ < 3.5 by incorporating 
the Forward Phig Calorimeter used during the HERAI running period^. Meaurements 

of cross sections as functions of Q^, xbj, E^t\ab' ^lab reasonably well-described by 
DGLAP-based calculations from Disent^ with large theoretical uncertainties at both low XBj 
and high rfl^^. Predictions from Lepto (DGLAP); AriadneHS! (CDM); and Cascade, with 
two sets from the J2003 unintegrated gluon PDF used, were also compared to the measure- 
ments. Overall, Ariadne provides the best description of the measured cross sections; Lepto 
consistently underestimates the cross sections, and Cascade fails to consistently reproduce the 
shapes of the distributions (see Fig. [3]). 
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Figure 3: ZEUS forward jets as a function of the kinematic variables xsj and compared to 
predictions from Ariadne, Lepto, and Cascade 



5 Trijet Production and Correlations 

Trijet cross sections and correlations were measured 
by the HI collaboration study of parton evo- 
lution at low XBj^^. Cross sections were measured 
as functions of xbj, jet pseudorapidity, scaled jet en- 
ergies, and correlations in the jet angles 9 and ip . 
The variable 6 is defined as the angle between the 
proton beam and the jet with the highest transverse 
energy, while ^ is defined as the angle between the 
plane defined by the proton beam and the highest 
Et jet, and the plane defined by the two jets with 
the highest Et- These measurements were made 
for three separate trijet samples: an inclusive trijet 
sample, and two trijet samples with one and two 
forward jets, respectively, with a forward jet having 
Cab < 20° and xjet = E^ncM/Epbeam > 0.035. 
For the inclusive trijet sample, NLOjet calcula- 
tions provide a reasonable description of the mea- 
sured cross secitons, but slightly undersestimate the 
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Figure 4: HI trijet cross sections as a function of 
xb] for events with two forward jets compared to 
NLOjet calculations at 0{al) and 0{al). 



measurements in the lowest bin of xbj • The agreement between the calculations and the mea- 
sured cross section in xbj is worse for the trijet sample containing two forward jets, with the 
most noticable disgreement observed at lowest xbj (see Fig. H]). The selection of two forward 
jets favors events with forward gluon emission unorderd in transverse momentum, which the 
calculations at 0{a^) do not predict entirely. Also seen in Fig. [J] is that the higher-order terms 
in the NLOjet calculations are important for forward jet emissions. The other cross sections 
for this sample are well-described by the calculations. 

Predictions from DjANGOH (CDM) and Rapgap LO MC models were also compared to the 
measured cross sections. The cross sections for the inclusive trijet sample are better described 
by CDM predictions, but both the CDM and Rapgap predictions are inconsistent for the jet 
correlation angles 6 and ip ; the Rapgap predictions fail to describe the 9 distributions, and 
the CDM predictions fail to describe the cross sections. 

6 Summary 

Parton dynamics at low xbj (10~^ < xbj < 10^^) have been investigated at HERA by the ZEUS 
and HI collaborations. DGLAP-based NLO calculations describe the measured cross sections 
and jet correlations reasonably well for the most part when higher-order terms in the calculations 
are properly taken into account. The calculations fail to describe the trijet cross section in xbj 
when the trijet sample contains two forward jets. Leading-order Monte Carlo models provide 
an inconsistent description of the measured cross sections. DGLAP-based LO MC models in 
general do not describe the cross sections; CDM models fail to describe Ai;^* and ip ; Cascade 
predictions are highly sensitive to the unintegrated gluon PDF used, and do not describe the 
data consistently. 
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